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Abstract 



Quantum error correcting codes enable the information contained in a quantum 
state to be protected from decoherence due to external perturbations. Applied to 
NMR, quantum coding docs not alter normal relaxation, but rather converts the state 
of a "data" spin into multiple quantum coherences involving additional ancilla spins. 
These multiple quantum coherences relax at differing rates, thus permitting the origi- 
nal state of the data to be approximately reconstructed by mixing them together in an 
appropriate fashion. This paper describes the operation of a simple, three-bit quan- 
tum code in the product operator formalism, and uses geometric algebra methods to 
obtain the error-corrected decay curve in the presence of arbitrary correlations in the 
external random fields. These predictions are confirmed in both the totally correlated 
and uncorrelated cases by liquid-state NMR experiments on ^"^C-labeled alanine, using 
gradient-diffusion methods to implement these idealized decoherence models. Quan- 
tum error correction in weakly polarized systems requires that the ancilla spins be 
prepared in a pseudo-pure state relative to the data spin, which entails a loss of sig- 
nal that exceeds any potential gain through error correction. Nevertheless, this study 
shows that quantum coding can be used to validate theoretical decoherence mecha- 
nisms, and to provide detailed information on correlations in the underlying NMR 
relaxation dynamics. 
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1 Introduction 



A quantum computer stores binary information in an array of two-state quantum systems, 
e.g. spin I nuclei. It performs logical operations on this information via unitary transforma- 
tions obtained by controlling the effective interactions among the systems. By operating on 
a coherent superposition over all combinations of states of the individual systems, it can, in 
effect, operate on all these combinations at once. This yields a degree of parallelism which 
grows exponentially with the size of the problem. Because such superpositions are extremely 
sensitive to decoherence, it was initially doubted that they could be maintained long enough 
to perform significant computations. Shor |l| and Steane were nonetheless able to devise 
error correcting codes to control decoherence in quantum computers. It has now been shown 
that arbitrarily long quantum computations can be carried out providing the error rate per 
operation is below some threshold 0, ^, |^, |^ . 

Even though a true quantum computer has yet to be built, small quantum computa- 
tions can be performed using standard liquid-state NMR spectroscopy to operate on the 
nuclear spins in ensembles of molecules in pseudo-pure states [|^, |[. These states may be 
characterized as having a density matrix with only a single nondegenerate eigenvalue, whose 
corresponding eigenvector transforms identically to the state vector of a true pure state. This 
approach has enabled all the basic operations of quantum computing to be demonstrated 

T^, including simple quantum algorithms [|TT|, |T2|, teleportation [r5|, and error correc- 



tion |15|. The preparation of pseudo-pure states from equilibrium states nevertheless entails 



a rapid loss in signal with increasing numbers of spins |Tq, O, 18], precluding the use of 



liquid-state NMR as a means of performing large-scale quantum computations. 

This paper builds on the authors' earlier implementation of a three-bit quantum error 



correcting code |15], and consists of two parts, one theoretical and the other experimental. 



In the first part, geometric algebra methods ||I9| are used to derive the error corrected decay 



curve assuming arbitrary correlations in the relaxation processes responsible. The second 
describes a complete NMR implementation of the code, and demonstrates that it performs as 
predicted by the theory using gradient diffusion methods to implement both totally correlated 
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and uncorrelated T2 relaxation. The theoretical (sections ^ — and experimental (sections |^ 
— parts of the paper may be read largely independently of each other. Section |2| provides 
an overview of the three-bit quantum error correcting code, while section ^ discusses the 
significance of the theory in the context of NMR. Together, the results of this paper show 
how quantum coding can be used not only to correct for decoherence, but also to design 
experiments to test specific hypotheses on the underlying relaxation mechanisms. 



2 The Quantum Error Correcting Code 

The theory of quantum error correcting codes, though only a few years old 0, @]' already 
a well-established and highly developed subject p, H, |, The "majority logic" code 
demonstrated in this paper protects information encoded in the joint state of three spins 
against arbitrary spin flip errors. Given a data spin (throughout this paper, the first) in 
some arbitrary state a|0) -|- (where |0) is the "down" (ground) state of the spin, |1) the 
"up" state, and |ap + = 1), two ancilla spins in their ground state are added to the 
system, obtaining (a|0) + /5|1))|00). A unitary transformation is then applied to the system 
consisting of two "c-NOT" (controUed-NOT) gates conditional on the data spin (see Fig. |l|). 
This step, called encoding, produces the state a|000) + 

Note that rotating a spin by vr in this state is the same as interchanging a pair of 
corresponding bits between |000) and |111), and results in an orthogonal state. Thus one can 
determine if a single spin flip has occured, and which spin it happened to, by a measurement 
that reveals which of the orthogonal subspaces a|100) -|- /?|011), q;|010) -|- /3|101), a|001) + 
/3|110) (if any) the system has moved into. If no error occured, this measurement leaves the 
state of the system intact, and otherwise one knows which spin needs to be rotated by tt 
in order to fix the error. This illustrates the basic idea behind all quantum error correcting 
codes: the correctable errors map the subspace of valid states onto mutually orthogonal 
subspaces, in each of which the unitary operation needed to fix the error is known. Because 
these entangled subspaces cannot be measured directly, the state must be decoded so that 
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the error syndrome is contained in the ancillae alone. This is done with the same unitary 
transformation as encoding, since it is its own inverse. 

In practice, the spins are subject to numerous small random perturbations from the 
environment arising, for example, from fluctuating external random fields, rather than single 
spin flips as above. Throughout this paper, these fields will be assumed to be about the x- 
axis in the rotating frame |^. To the extent that these perturbations affect multiple spins, 
the data spin will lose coherence in the yz-plane despite error correction. It nevertheless 
follows from the linearity of quantum mechanics that the code will still prevent decoherence 
to first order in time, as will now be shown. 

The interaction of such fields with the y and z-components of the spins has the Hamilto- 
nian 

Hx = I'Bimi + ^'Bm^i + I'B'^mi , (1) 

where denotes the x-component of the fluctuating magnetic field at each spin, is its 
gyromagnetic ratio, and = |<t^ are the usual angular momentum operators in units of h 
{k = 1, 2, 3). The corresponding propagator has the form exp(— + + where 
L denotes the imaginary unit and x'^ = 1^ Jq B^(T)dT. In the limit of short decoherence times 
and hence small accumulated phase errors x'^; ^ first-order expansion of this propagator is 
adequate, namely: 

exp {-cixX + XX + XX)) ^ 1 - ^XX - ^X'll - ^XX (2) 
Thus to first order this random propagator maps the encoded state to 

(a|000) - fx'("|100) + /5|011)) 

- fx'(a|010) + /5|101)) (3) 

- fx'(«|001) + /5|110)) . 
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The decoding operation maps this to 



(a|0)+/5|l))|00) - fxH«|l) + /?|0))|ll) 

- L^^(a\0) + pmiO) (4) 

- fx'(a|0) + /?|l))|01) . 

Finally, a Toffoli gate (also called a "c^-NOT") corrects the error, by flipping the data spin in 
the term in which both ancillae are "down" This returns the data spin to its original 

state a|0) + in that term, and otherwise leaves it alone. 

Thus, in addition to single spin flips, the code corrects to first order for random, time- 
dependent rotations about the x-axis as claimed. By inserting a 7r/2 y-rotation of all the 
spins after encoding and its inverse before decoding, this code can be transformed to one that 



protects against random z-rotations, which give rise to adiabatic T2 relaxation in NMR [[Tq] . 
Because the code does not protect the ancillae, they must be "discarded" after correction, 
which corresponds to taking the partial trace over them. In NMR spectroscopy, this operation 
is performed by decoupling the ancillae from the data during acquisition. Note also that, 
since strong measurements are not available in ensembles, it is not possible to reset and reuse 
the ancillae for repeated error correction (see section ^ for further discussion of this issue). 

The following three sections will present a detailed theoretical analysis of the code in the 
product operator notation of NMR [ ^0] , p5| , using geometric algebra methods []T^ to compute 



the results of the encoding, decoherence, decoding and correction operations. Readers who 
are interested primarily in the experimental procedures and results should proceed directly 
to the main theoretical result, Eq. (^5[), and continue reading from there. 



3 Geometric Algebra Analysis of Encoding 



In NMR spectroscopy, one does not observe individual spin systems, but rather averages 
over macroscopic ensembles of such systems. In the semi-classical treatment of liquid-state 
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NMR relaxation processes used here PT| , two distinct types of averages must be taken. The 



first is an average over the initial state of a typical molecule, which results, for example, 
in the equilibrium state + + I^. The second is an average over the environment of a 
typical molecule, or more precisely, the external random fields which act upon it. Provided 
the random fields are uncorrelated with the initial state, these two averages can be taken 
separately. In order to focus on how the results of error correction are averaged over the 
environment, in the following three sections it will be assumed that the spin system is initially 
in a pure state; the effect of averaging over an ensemble of spin systems in different pure 
states is considered in an Appendix. 

Even though an initial pure state is assumed, the random fields from the environment will 
cause the state to decohere into a general mixed state, which must be described by a density 
matrix. The physical meaning of a density matrix is most clearly expressed using the product 
operator formalism, since it is based on a complete system of observables |]22|, ^ . In the 
algebra generated by the product operators, spinors (pure states) correspond to idempotents, 
namely elements that are equal to their squares. In particular, spinors of the form 1 5^5^(5^) 
(5^ = 0, 1; = 1, 2, 3) correspond to products of three idempotents, namely: 

E^E|Ei = 1(1 + (-l)^'2li)i(l + (-l)^' 212)1(1 + (-l)^'2l3) ^ |5W)(5W| (5) 

(Note that throughout this paper the scalar identity 1 is identified with the 3-spin identity 
operator Ig.) Letting a, P be the complex conjugates of G C, and 3?, denote the real 
and imaginary parts, the density matrix of the data spin #1 in a general superposition state 
a\0)+p\l) is: 

pi ^ (a|0)+/?|l))(a(0|+/3(l|) 

= i(|«r + |/?r)(|0)(0| + |l)(l|) + J?(«/3)(|l)(0| + |0)(l|) + 

.9(a/3)(|l)(0| - |0)(1|) + i(l«P - m{\0){0\ - |1)(1|) (6) 

= + i/^n + + ^i^p) 2i; + (i^r - ii 
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Alternatively, an arbitrary state can be obtained by rotation of the |0)(0| = state, which 
may be specified in terms of polar angles as 

= e-'-^^^'(cos(|) - 2tsin(f)Ii)E^cos(|) + 2z, sin(f )Ii)e^'^i^' (7) 
= (cos(l)e-^^i^' - 2tsm{l)lle"f'^i)E\{e''t'^"^ cos(f) + e-^'^i^'2t sin(|)Ii) 
= (cos(f)e-^^/2 _ 2i,sin(|)Iie^<^/2)E^(e^<^/2cos(|) + e-^<^/22tsin(f )Ii) , 

where the last line was obtained by using the relation I^E^^ = ±|E]i_ to absorb the ll 
operators in the exponentials. A comparison of these two expressions for p\ shows the 
coefficients in the superposition may be interpreted geometrically as 

a= cos(|)e-^*/2 , (3= - 1 sin(f ) e^*/^ . (8) 



This correspondence between states and operators, which is the basis of the product 
operator formalism, has been used together with reformulations of multiparticle quantum 



mechanics in the language of geometric algebra |2^, ^ |2^, Q to obtain algebraic expres- 
sions and geometric interpretations for the basic operations of quantum computing |19|. Of 
particular interest is the fact that the propagators of c-NOT gates can also be written in 
terms of idempotents. Thus the particular c-NOT gate that flips the second spin when the 
flrst is "up" has the form: 

S^li = exp(t7rEi(l-2l2)/2) = Ei exp(t7r(l-2I^)/2) + (l-Ei) (9) 
= 2I^Ei + E^ = 1-Ei(l-2I^) 

Using the relations 41^12 = 1, i2e2 = l2 and E^EL = E^E^ = 0, S^l^ may be shown to 
exchange E^ with E^ in those products containing Ei; for example: 

g2|i(g^E^^)g2|i _ (2l2E^)(EiE^)(Ei2I^) +E^(EiE^)E^ 

= E^2I^i(l + 2I,2)2I^) (10) 
= E^i(l -SI^I^I^) = E^E^^ 



In a similar fashion, it can be shown that the Toffoh gate, which flips the first spin only if 
the other two are "up" , has the form: 

T^|23 = exp(t7r(l - 2I;^)E^E^) (11) 
= 2I^E^E^ + (1 - E^E^ ) = 1 - (1 - 2I^)E2 E^ 



These are all the gates that are needed to implement the error correction procedure studied 
in this paper, whose operation will now be analyzed using geometric algebra. 

The first step of the procedure is to encode the state p\ of the first spin by correlating 
it with the ancillae in the state E^Ei]_. This is done by applying the c-NOT's S^'^ and S^'^, 
as shown in Fig. |l|. These c-NOT's commute, and their product is 

S2|is3|i _ (2l2Ei + E^)(2I^Ei + E^) (12) 

= 4I^i3eL + E^ = S^^li . 

It follows that 

S23|i2iis23|i ^ (4I2i3eL + E^)(2I^)(4I2i3e^ + E^) (13) 
= (2li)(4l2lX + Ei)(4I^l3Ei + E^) 
= (2Ii^)(4I^I^(EUE^)) = 8I^I^I^ 

Since S^'^'^S^'^'^ = 1 and S^'^'^ commutes with Ej'^E^E'^, the density matrix pg of the encoded 
state, which is obtained by applying S^^'^ to the initial density matrix = p^E^E"^, is 
given by 

= S2=^|i(a + f32li) S23|iS23|i(E^E^E^)S2=^|iS2=^|i(a + P2ll) S^^l^ (14) 
= (a + /38I^I^I^)(E^E^E^)(a + /^SI^I^I^) . 

The product of all three Ix operators in this expression refiects the entanglement of the three 
spins in the encoded state. 
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4 Averaging Over the Environment 



The next step is to apply a random rotation about the x-axis to the encoded state, and then 
to compute the average of the result. Fortunately, the x-rotation commutes with the outer 
factors of pg, so that it can be applied directly to the inner factor: 

Since the outer factors are constant, they can be taken out of the ensemble average. This 
reduces the analysis of the decoherence process for arbitrary a and (3 to the special case in 
which a = 1 and (3 = 0. 

To facilitate this process, rewrite the inner transformation as 

= 1(^1 + e-2^^'i^'2Ii) (l + e-2^^'^'2l2) (l + 6-2^x^^.^213) . 
Expanding this product yields the sum of all possible terms of the form 

Rsu^s^x', X\ X') ((1 - S') + 6'2ll) ((1 - 6') + 6^211) ((1 - 6') + 6'2ll) (17) 

where 

Rs^sMx\x',x') - e--(^^-^^^^-^^^-^^^^+^^^^^^^) (18) 

for 6\6^,6'^ e {0,1}. 

It follows that the average state can be obtained by applying the averages of the propaga- 
tors R5i5253(x^, X^) to the corresponding terms of the encoded state pg. Assuming that 
the spectral density of the random fields B^(t) is well-approximated by a delta-function (i.e., 
that the correlation time is small compared to the inverse square root of the variance in the 
frequencies ■jBy.(t); see for details), the joint probability density function of the random 
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angles x\x iX will be a multivariate Gaussian whose covariance matrix grows linearly in 



time. Letting x = [x^? X^? X^V = [c'''t] = Ct be the covariance matrix, this 

can be written as 

P{X) = ((27r)3 det(Ct))"^e-5t>^^<^"''^ . (19) 

At this point it turns out to be easier to rotate 'R.s'<-5'2S3{x^, X^ ■, X^) to the z-axis: 

KhhAx\ x\ x') = e^f (^H^^^?)R,.,.,3(x\ X^ x')e-i(^^^+^?+^?) 

_ ^-2,{s^x'ii+s\^i^.+s\'i!) ^20) 

This is because the idempotents E± "absorb" the Iz operators in the exponent, so that it 
can be expanded into a sum of terms involving only scalar exponentials, 

R'.l...3(x^x^x') = R's^s^s^.{x\x\x'){K + ^-)iK + ^-)iK + ^'^ (21) 

where each term has the form 

R's.s^M\x^x')K^^%n-^ = e-(^^>^^+^VH-.3x^)Ei,E^,E^3 (22) 

with e'' E { — 1, 0, +1}, e'^ = S'^e'' and Eq = 1 for A; = 1, 2, 3. The average of each such scalar 
exponential is 



/oo 
P(x)e-^'^-^dx = e-t^"^^ , (23) 
-oo 



where e = [e^,e^,e^]^. Since this is independent of the overall sign of e, for e ^ the 
ensemble average of each R^i^a^slxS X^; X^) is a sum of pairs of terms of the form 

e-l-^Ce^gi^E^.E^a + Ei^iE^^.E^^a) . (24) 

In the case that = = = 0, Eq. (^3]) yields 1 as expected, while in the case that 
= ±1 while the other two e* = (1 < i 7^ j < 3), the two terms in Eq. (P^ ) are 



-2^Ji = e-5^"(Ei + Ei) = e-*'^"/^ (25) 
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In the case that e' = while and e'^ are nonzero {{i,j,k} = {1,2,3}), Eq. (P^ ) has the 
form 



^_|(,..+e-+2.^.'=c^'=)^E,.E.^ + E_,.E_,.) (26) 

where 

(E,.E,.+E_,.E_,.) = i(l + 4e^-e^F,I,^) ^ E% (27) 

and e!^ = 1(1 ± 4ip^) is idempotent. Since there are two such pairs of terms, one with 
= and the other with = —e^, the complete result is 

g-2.(xii^+x'=iz^) = e"t(^''+^'') (^e~*'=''Ej + e*'''Ei'') 

= e-5(^''+^"') (cosh(t d^) - sinh(t <^^) Al{l^) (28) 

= g-|(c^'^+c'''=+8cJ*I^'lJ) 

In a similar fashion, it can be shown that the average propagator for random rotations acting 
on the 81^1^1^ term is 



g-2t(xlll+x2l2+x3lt) = g-|(c"+c22+c33+8cl2lll2+8cl3lllf+8c23l?l|) j^29) 

These averages are easily rotated back to the average of the original propagator (x^, x^? X') 

simply by replacing by Ix throughout. All the terms obtained by the above expansions 
are collected in the next section. 



5 Decoding and Error Correction 

In order to write the ensemble average of Eq. (0) in a compact form, define the time- 
dependent (nonunitary) operators 

= F^(t) = e-*""/2 (30) 

(acting by scalar multiplication), and 

pifc = F{\t) = e-*^''4«i' (31) 
10 



(acting by left-multiplication). Then the density matrix after decoherence becomes 

Pc ^ (a + PSIIIX) Tc [1(1 + FH\^:){\ + F22I,2)(1 + FH\% (a + pSiXll) , (32) 
where is a linear operator- valued function defined on the products of operators by 

^c[l] = 1, 

^c[2Fj = 2Ii (l<j<3), 
J^cWz] = Ff4Fj,^ (l<i<A;<3) (33) 

= g-tci24lil2-tci34lil|-tc234l2l3gjlj2j3 _ 

Next, the decoding operation is performed by applying the same two c-NOT's used in 
encoding, i.e. S^^'^. On inserting appropriate factors of unity, the resulting density matrix 
can be written as 

Pd = (a + (32ll)J^o [| S23|i(l + F'2Il)(l + F^2Il)(l + F'2Il) S^^li] (a + P2ll) (34) 
= (a + P2Jl)J^o [I (1 + F'21l){l + FHlllDil + FHIII^] {a + /32I^) , 
where the linear operator- valued function J^d is is defined on the products of operators by 

This translates to 

J^d[1] = 1, 
To[2ll] = 2ll, 

J^o[2ll] = Fl,Hll = e-''"'^'i'^2ll, 
11 



= 41^1^ 



(36) 



To[8lim ^ F^2F^3f23 8J1J2J3 



^ g-tel24lllj-tei341il2-ic234l2l3 gjlj2j3 



since 



g-tci2 41^13 _^^^_2); 

g-tc234l2l3 ^j^2,k = 3). 



Finally, the error correcting Toffoli gate T^l^^ = 2I^E^ + (1 — E^E^ ) is applied. Since 
'pi|23 commutes with the outer factors a + /?2I^ and its conjugate, one could define a new 
function and work out how it operates on the transformed states as was done with JF^ above. 
Rather than simply permuting product operators as S^^li does, however, T^'^^ maps each 
to a linear combination of four products, and in order to evaluate the partial trace over the 
ancillae it is necessary to fully expand the result. Fortunately, this laborious task can be 
avoided by noting that for any product operator expression X, the partial trace over the 
ancillae 4(X)23 jg g^j^g ^g 4(£:[x])23 

= 4 (E^E^ X E^E^ + E^Ei X E^E^ + E^E^ X E^E^ + E^E^ X E^E^ (38) 

(This follows, for example, from the product operator characterization of the partial trace 
given in Ref. Eq. (40).) In addition, the Toffoli commutes with all of the idempotents 
Ej. (j = 2,3), so that S can be applied before applying the Toffoli gate. As shall be 
seen momentarily, this dramatically simplies the intermediate results, so that the linear 
combinations created on applying the final Toffoli are easily dealt with. 
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Because the idempotents commute with for all 1 < < j < 3, we can apply S 
directly to the coefficients F-^'^ in Eq. (|36|) , e.g. 

^[Fd] = E22E^3(cosh(tci2)-4I^I^sinh(te^2))E22E33 (39) 

= (cosh(tc-) (E^.E^a)^ - 41;^!^ sinh(tc-) (E^.E^.^a) (E^.E^a)) 

e2,e3e{±l} 

= COsh(tel2) Y E'^E'^ = COsh(tcl2) ^ pl2 

e2,e3e{±l} 

In a similar fashion, it can be shown that 

F^^ = ^[F^=^] = cosh(te^3) , = S[Y^^] = cosh(te23) , (40) 

and = 

^jpl^2pl3p23] = COsh(tel2) COsh(tcl3) COsh(tc2=^) - Sinh(tc^2) gi^l^(^^13) gi^^(^^23^ _ (4^^ 

Thus applying the Toffoli to the projection of each term in Eq. (^) yields 





= 1, 


Tl|23f[J-D[2I,l]]Tl|23 


= + 2l\ll + 21,^13 - Al\llll , 






ti|23£:[^d[2i3]]ti|23 


= 21^ , 


T1|23^[J-d[4I^I,2]]T1|23 


= + 21^1^ - 21,^1^ + A\\llll , 


T1|23£[J-d[4I^I3]]T1|23 


= 1\- 2l\ll + 2I,il3 + Al\llll , 




= 41^1^ , 


T1|23^[^d[8I,1i2i3]]t1|23 


= (-1^ + 21^1^ + 21^1^ + 



The partial trace over the ancilla simply drops all terms from the above that contain factors 
of or I^, which eliminates all lines from the above equation containing a doubly-indexed 
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time-dependent coefficient F^^. Multiplying each of the surviving terms by the appropriate 
scalar exponential = exp{—tc'^/2) and combining finally yields 

Pi = {Per - {{a + P2ll)T'^''p^T'na + P2ll)y' 

= {a + P211)(^IT^\^''S[J^d[1 + F'2II + FHiIiI + FHiIiI (43) 

+ F'F^F^ Fi23 8jij2j3 ] ]Ti|23^'' (« + ^ 21^) 

= {a + 13 2ll) 1{1 + {F^ + F^ + F^- F^F^F^ F^^^)ll) (a + /3 21^) 

= 1 + ^{a(3)2ll + (3(a/3)2i; + {\a\' - \(3\')ll) e(t) , 

= 1 ( e-*-"/2 + e-''^""/^ + e-*'^''/^ ) _ i e-tc"/2 ^-tc^V2 ^-tc^V2 ^ ^44) 



where 



This equation describes how the error-corrected state of the data spin #1 decays in the 
presence of decoherence due to external random fields about the x-axis, where tc^^, td'^ are 
the variances and covariances among the phases of the three spins, as described in section ^ 
(Eq. dg)). 

It is readily verified that the derivative 0(t) vanishes at t = regardless of the covariances, 
as expected from the analysis given in section ^ Two special cases are of particular interest 
in the following. The first occurs when the random phases ^ire identically distributed and 
uncorrelated, so that c^^ = c^^ = c^^ = 2/r and c^^ = c"^ = c^^ = 0; then the above simplifies 
to: 

e(t) = i ( 3e-*/^ - e-=^*/" ) (45) 

The second occurs when these random variables are totally correlated, so that c^'' = 2/r for 
l<3 <k<3; then: 

e(t) = I (9e-*/"-e-^*/") (46) 
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6 Discussion of Theoretical Results 



Quantum error correcting codes were designed to prevent decoherence from occuring within 
a subsystem, at the expense of increasing decoherence in a quantum environment with which 
it interacts in a controllable fashion. In order to devise a quantum error correcting code, the 
mechanism of decoherence must be precisely known. Since a code can only correct for decay 
due to decoherence to first order (in the absence of concatenation [^]), decoherence can only 
be prevented over long time intervals if the error correction procedure can be performed 
much more rapidly than decoherence occurs. To do this, it is necessary to be able to rapidly 
reset the ancillae to their ground states after each correction so they can be reused, or else 
to have a large supply of ancillae in their ground states that can rapidly replace the used 
ones. Neither of these conditions is true in liquid-state NMR, since Ti > T2 and the number 
of spins whose interactions can be precisely controlled is limited. 

Of course, pure states are not available at all in liquid-state NMR spectroscopy, as was 
assumed above for the ancillae. The isomorphism which exists between the dynamics of 
pure and pseudo-pure states ensures that error correction will nevertheless work, so long 
as the ancillae are in a pseudo-pure ground state relative to the data spin. As discussed 
in detail elsewhere [|^, the preparation of the ancillae in a pseudo-pure state entails 
a loss of polarization that exceeds any potential gain through error correction, so that the 
net signal-to-noise in the NMR spectrum is actually lowered by the use of error correction. 
It would therefore be advantageous if a broader class of mixed states for the ancillae were 
also suitable for error correction, but arguments are given in an Appendix which imply that 
this is probably not true, and certainly not for any diagonal (z-polarized) mixed state of the 
ancillae. 

In NMR, where decoherence occurs via a Raman process involving both translational 
and rotational molecular motions, measurements of the covariances in the fluctuating fields 
are important indicators of the nature of these motions. Thus, even though error correction 
cannot be used to sharpen the lines or improve signal-to-noise in NMR spectroscopy, the idea 
of coherently mixing together single and multiple quantum coherences so as to cancel the 
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first-order decay, under certain specific assumptions regarding the underlying mechanisms 
responsible for the relaxation and the correlations among them, clearly has the potential of 
being useful in NMR studies of molecular motions. The higher moments of the initial decay 
curve should also provide further confirmation of the assumptions. For example, the second 
derivative e(t) at t = is 9(0) 

= -i(2(ci2)2 + 2(ci3)2^2(c23)2 + ciic22 + c"c33 + c2V3) <0, (47) 
while the third derivative is G(0) 

= ^ (3(c")2(c22 + c33) + 3(c22)2(c" + c33) + ^{c'Yic'^ + c'') + 6{c''c''^c'') 

(48) 

+ 12((cl2)2^ (^13)2 ^(^23)2)(^11^ ^22 ^^33) ^48(^12^13^23)) ^ 

which we expect will be generally positive. 

Decoherence is of central importance to quantum mechanics, since it is the process by 
which classical statistical mechanics emerges from the underlying deterministic evolution of 
wave functions |3^, 0. It is also a process which has proven extremely challenging to study 
in its full generality |3^, Although again somewhat limited by its ensemble nature, 
NMR is both theoretically and experimentally a very convenient model system in which to 
study decoherence. Because quantum codes only protect against decoherence by specific 
mechanisms, they can be used to design NMR experiments to test the validity of theoretical 
models. This will now be illustrated by describing liquid-state NMR experiments which 
validate the forgoing theoretical results. 

7 Gradient Implementation of Decoherence 

Loss of phase coherence is, of course, a natural relaxation process in NMR spectroscopy. 
The longitudinal fiuctuating fields that induce this decoherence, however, come from many 
different sources, both intramolecular and intermolecular, and involving both spin and non- 
spin degrees of freedom. Superimposed on these longitudinal fluctuations are transverse 
fluctuations, which involve an exchange of energy with the environment. The three- bit 
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code majority logic analyzed above can deal with fluctuations about only one transverse 
axis; by including 7r/2 rotation and its inverse at the end of the encoding and beginning of 
the decoding, respectively, it can deal instead with longitudinal fluctuations. (Fluctuations 
about all three axes can only be corrected using at least four ancillae 0, ||, |], In 
addition, intramolecular dipole-dipole interactions are an important source of decoherence 
in liquid-state NMR, and these will also not be corrected by the three-bit code. Thus the 
forgoing theoretical results can only be rigorously demonstrated by experiments in which the 
dominant relaxation mechanism is a known, artificially induced fluctuation about a single 
transverse axis. 

This is done here using molecular diffusion to randomize the spatial variations in the 
phase created by a pulsed field gradient. A field gradient along the z-axis causes the spins 
to precess at different rates, depending on their z-coordinates. This winds the transverse 
magnetization due to each coherence into a spiral, whose pitch decreases linearly with the 
length of the gradient pulse. In a liquid sample, the decay of a tightly wound spiral is due 
almost entirely to diffusion along its axis. Because the change in phase due to diffusion is 
exactly the same for every spin in a molecule, this implements the totally correlated error 
model. 

The quantitative analysis of this process is best done using the fc-space formalism [^j . 
The Fourier transform of the distribution of transverse magnetization v{z) along the z-axis 
will be denoted here by T{k). The spiral produced by the gradient is described by the 
product v{z) exp{tnijkoz), where riij is the order of the coherence in question, whose Fourier 
transform is T{k)6{nij{k — ko)). Diffusion operates on the magnetization distribution by 
convolution with a Gaussian whose variance is cr^ = Dt, where D is the diffusion coefficient 
3^ , i.e. v{z) expiiUijkQz) -kexp^—z"^/ (2cr))/-\/27ro^. The Fourier transform is thus multiplied 



by the Gaussian exp(— (A; — ko) nqjD/2), so that an inverse gradient pulse produces an echo 
which decays with the square of the coherence order, as expected from a totally correlated 
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dephasing process. These relations may be summarized as follows: 
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/V2 



(49) 



T{k) g^^^'""^ T{k)6{n,,{k - k,)) diffusion ^^^^ ^ ^^^^^.^^ _ k,))e-^l''^'/' 



Implementation of the uncorrelated error model was done by applying three successive 
gradient- diffusion sequences. In each sequence, two of the spins were refocussed during a 
gradient-pulse by selective RF vr-pulses, while maintaining the phase ramp on the third spin 
(which was different in each of the three sequences). This was followed by a time interval 
to allow diffusion to randomize the phase of the third spin, after which a second gradient 
and further selective vr-pulses were used to refocus all three spins (save for the coherence of 
the third lost to diffusion). This procedure randomized the phase of each spin equally and 
independently during separate diffusion intervals, as required by the uncorrelated model. 
The complex RF pulse and gradient sequences which performed this task are described in 
the following section. 



8 Experimental Implementation of Error Correction 

This section provides a detailed description of the RF and gradient pulse sequences used to 
implement the three-bit quantum error correcting code (see Fig. |l]) as well as the correlated 
and uncorrelated decoherence models. The challenge here lies in obtaining precisely the 
desired effective Hamiltonian at each step, and in putting these steps together without side 
effects. In particular, the finite duration of all experimentally realizable RF and gradient 
pulses allow the system to evolve and so pick up wanted changes in the relative phases of 
its states. Therefore these unwanted evolutions must be refocussed, at the cost of increased 
complexity in the implementations. Additionally, the system is subjected to incoherent er- 
rors due to pulse imperfections and field inhomogeneity. These could be compensated for by 
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phase cycling and other averaging techniques but this is undesirable in the present 

context since it complicates the interpretation of the experiments. Of course, such difficulties 
are not unique to NMR, but are expected to varying degrees in any quantum information pro- 
cessor. The intrinsically long decoherence times, innate averaging over ensembles and superb 
coherent control available in modern NMR spectroscopy is what makes these experiments 
possible today. 

The three-bit quantum error correcting code was realized with a sample of ^^C-labeled 
alanine (NH+-C"H(C'^H3)-C02 ) in D2O solution at room temperature. The measurements 
were carried out on a Bruker AMX400 spectrometer (9.6 T) equipped with a 5 mm probe 
tuned to the ^'^C and frequencies of 100.61 and 400.13 MHz, respectively. This probe 
was equipped with a z-gradient coils capable of generating field gradients of 60 G / cm. 
With decoupling of alanine's protons, this system exhibits a weakly-coupled carbon NMR 
spectrum. The internal Hamiltonian in the rotating frame is 

Hi,, = o^C'lf +,;C^lf + 27rjC'C"lflf + 27rjC'C''lflf +27rjC"C^lflf (50) 
where (with the transmitter on-resonance with the C") 

uj^'/{2ti) = 12580Hz, cj^7(27r) = OHz, lj^^/{27i) = -3443Hz, (51) 
JC'C" = 54.2Hz, JC'C" = 1.2Hz, J^"^' = 35.1Hz . 

The was chosen as the data spin 7^1, and exhibits a well-resolved quartet in the spectrum. 
This has the advantage that none of the gates required for the error correction procedure 
needed to use the small 1.2 Hz coupling directly (although it was necessary to refocus this 
coupling). The carbonyl C (spin #2) and (spin #3) were used for the two ancilla spins. 

To implement a universal set of quantum logic gates in a spin system it is sufficient 
that all the spins be connected by coupling pathways and sufficiently well-resolved to enable 
arbitary rotations to be applied to any subset of the spins. Here selective excitations were 
implemented using phase-modulated Gaussian pulses |^. These shaped pulses consisted 



of 384 complex points with a total duration of 1.5 msec, with an excitation profile that 
consisted of Gaussians centered on the frequencies of the (one or more) spins of interest with 
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a standard deviation of 500 Hz. The transmitter was placed on the data spin, enabhng it to 
be controlled by soft rectangular pulses. 

It is convenient to define five specific pulse sequences, which were designed to yield a 
simple, well-defined effective Hamiltonian and served as the "modules" (building-blocks) for 
the overall pulse sequences. In the following list, all the time intervals are calculated from 
the midpoints of their surrounding shaped pulses. 



Identityit): This is a "time-suspension" sequence which is designed to refocus all components 
of the Hamiltonian by the symmetric application of vr-pulses . Identity sequences correct 
the phase errors due to the finite duration of RF and gradient pulses, and are used whenever 
two noncommuting operations must be applied sequentially. The implementation used in 
these experiments is given by the pulse sequence (in left-to-right temporal order) 

A - vr - -TT - -TT , (52) 



TT — TT — — TT 



where {tt)^^ denotes a vr-rotation of spins about the x-axis, and similarly for (vr)^™ and 
(tt)^™. These pulses have a duration of A, so that the net propagator can be written as e.g. 

with similar expressions for the other pulses. Thus if one places such an identity sequence 
between two noncommuting pulses, the evolution during the A/2 duration of the last half of 
the first pulse and the first half of the last pulse (as in Eq. ([53|)) cancels with the evolution 
during the intervening identity sequence. Note that due to the symmetry of this sequence all 
the CH couplings are refocussed, so that additional decoupling of the protons is not needed. 

Jdelay{k,i,t): This module yields the effective Hamiltonian 27rJ'^'l^I^ where k,l = 1,2,3 
are spin indices. The parameter t determines the net acquired phase = 2'kJ^H\ if = 1, 
£ = 2, for example, a phase shift of is obtained from t = (p/ (2 J^^) = ip/{2 x 54.2) sec. The 
complete sequence is 
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where [A] is again a delay of duration equal to that of the shaped pulses (1.5 msec). The 
finite duration of the pulses sets a lower limit for the phase evolution of > 4 ■ 2 ■ J^^A 
(= 0.2l7r for J^^ and O.lSvr for J^^). 

Jdelaylnv{k, i, t): This is a short version of J delay with one modification: after its application 
all three spins are inverted by vr relative to the result of Jdelay. Successive evolutions under 
27r J^^I^I^ and 27rJ^^I^I^ were often needed in these experiments, and the same result is 
obtained from using two successive Jdelaylnv modules as from two successive Jdelay modules. 
The former, however, saves three shaped pulses, as may be seen from its pulse sequence: 

(0-(-)"-a)-(-)"-(i)-(-)"-a) ™ 

TC-Decohere{g, 6, t, T): This module was used to induce decoherence under totally correlated 
external random fields. It includes two gradients of equal strength g and duration 6 but 
opposite polarity, embedded in a time-suspension sequence of length T. The gradients are 
separated by a period t during which diffusion takes place. The effects of decoherence due to 
T2 relaxation during this module were reduced to a constant factor by keeping the total time 
required fixed at T while varying t. Letting [z— grad(5', 6)] denote a gradient pulse along the 
z-axis of strength g and duration 6, the pulse sequence used for this module was 

- (-g-dto.*)) - (5 - 5) - (^) -{'2- 1)- (-s-d(-«,«)) 

This pulse sequence is shown in diagrammatic form in Fig. |^(a). 

UC-Decohere{g,6,r2,ri): This module was used to induce decoherence under uncorrelated 
external random fields (see above). It uses a combination of refocusing and gradient pulses 
of absolute strength g and duration 6. These gradient pulses were strung together into one 
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of three sequences, all of the form 

G': (^z-grad(^?i, 5)^-0 - (^z-grad(^?^, 5)^ - - (57) 

z-grad(^^, 6)^ — ^ — 0grad(c/^, S)^ — (^n 
where the polarities in the sequence for dephasing the z-th spin were given by 

9^ ■ [+9, -9, +9, -9] 

9^ ■■ [+9, -9, -9, +9] (58) 

9^ ■■ [+9, +9, -9, -9] ■ 

Each of these gradient sequences, in turn, was embedded in a fixed number of repetitions of 
a sequence of refocusing pulses, of the form 

^)-0 -0-0 ' (59) 

where (5) indicates a time delay. Note that two consecutive repetitions of this sequence 
refocuses both the chemical shift and coupling evolution, and that the sequences above 
with g = comprise two such repetitions and so do nothing, as desired. A given number R of 
repetitions of this pulse sequence will be denoted by (R), while the gradient sequence which 
refocuses the effects of the gradient sequence (G*) in a given experiment will be denoted by 
(G*). Depending on whether the intervening number of refocusing sequences ( |59D is even or 
odd, this will either be (G*) or the same sequence with all its polarities reversed. Thus the 
overall sequence used to independently decohere all three spins and thereby implement the 
uncorrelated decoherence model is given by 

r2 - n] (g'] (n) (g'] (r, - n] 0^ (60) 

=(3 



n - G' - r,-n - - n - G' 



where < ri < r2 are integers and (ri) or (r2 — ri) denote the enclosed number of repetitions 
of (p9D. This pulse sequence is shown in diagrammatic form in Fig. @(b). Observe the total 
length of each experiment is constant, which keeps the effect of intrinsic T2 relaxation to a 
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constant overall factor. The signal loss due to RF inhomogeneity during the application of 
the long train of refocusing pulses is reduced by a CPMG type phase modulation of length 
eight [0. These phase changes are carried out during a single scan of each experiment, and 
no phase cycling is performed across scans. 

In the TC-Decohere module, the gradient strength and duration were g = 35.7 Gauss/cm 
and 6 = 2.5 msec, respectively, while the total time was kept fixed at T = 64.5 msec, with 
a diffusion time increment of t = 4 msec. In the UC-Decohere module, the gradient strength 
used was \g\ = 12.2 Gauss/cm, with 6 = 2.078 msec, which yields an increment of 7.156 
msec, over r2 = 21 time points. 

An account of how to design pulse sequences, composed of these modules, for any desired 
effective Hamiltonian may be found in our earlier work [|l^, For the sake of complete- 
ness, however, the basic ideas are repeated for the case of the error-correcting Toffoli gate. 
This gate may be written in exponential form and expanded into a product of commuting 
propagators, as follows: 

r[,i|23 ^ 2I;^E^E3 + (1 - E^El) = e^'^(i/2-iiEiE3 )e2_e3_ ^g^^ 

^ gt^r/S g-«.7r/4ll g-t.7r/4l2 ^-Ln/4P^ g'.'r/2lll2 ^Ln/2I^^P^ ^'■^fiP^l^, ^,-''^1^^^! 

The propagators in this expression could be further expanded into products of propagators 
directly implementable using the above modules. In these experiments, however, the Toffoli 
gate is followed immediately by a partial trace over the ancillae, i.e. by observing the data 
spin while decoupling the ancillae. Hence all propagators that operate only on the ancillae 
can be eliminated with no effect on the final result. Since the net phase is also unobservable, 
this leaves only the product of the propagators 

g-67r/4li g*.7r/2lll2 g67r/2lll3 g-z-Trllj^lf _ j^g2) 

This in turn can be implemented by the sequence of propagators of the effective Hamiltonians 
(in reverse temporal order): 

gt7r/4li gt7r/2li ^lttIII^ ^-lttI] ^l-k/211iI ^lttI^I^ gi-7r/2li ^l77t/211I^ ^-l-kI] g*.7r/2lllf ^-ltt/21] ^^gg^j 
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These effective Hamiltonians can all be obtained by combining the above modules with RF 
pulses. 

The state I^E^E^ was prepared first, after which the other states I^E^E^ and IyE^Eij_ 
used in the experiments were obtained by rotating spin ^^1. Starting with the three-spin 
equilibrium state + + 1^, a pair of c-NOT's was applied to the ancillae conditional on 
the data spin: 

g2|is3|i (ji + j2 ^ s^l^S^li = + 2ll + 2ll) (64) 
This was then subjected to the pulse sequence (given in temporal order) 

(f )' - (47-) - (f)i/4 - (><-grad) - (f - (^) - (f)t/, - (y-grad) , (65) 

where [1/(27^*^)] is an evolution under the effective Hamiltonian 27rJ^'^I^I^ a time 1/(2 J^'^) 
[k = 2, 3), i.e. Jdelay{l, k, t), (x— grad) is a magnetic field gradient of dB^/dx, and similarly 
for (y— grad). This yields 3I^E^Eij_. Note that this method of preparation reduces the signal 
of the data spin by 25% (as opposed to the 50% expected from the equilibrium populations 
IITsI ); the transformation in Eq. ( |65l) also has the interesting property of being a projection, 
in that applying it a second time to the resulting state does not change the state. 

The final operation required is to observe the partial trace over the ancillae at the end. 
Because only single quantum coherences of the form and ly are observable, to implement 
the partial trace during acquisition it is only necessary to prevent antiphase magnetization 
with respect to the ancillae from evolving into observable single quantum terms. This was 
achieved by interspersing the sampling of the signal with vr-pulses on the ancillae, thereby 
repeatedly refocusing these terms. This RF irradiation caused a fixed upfield shift of ~ 12 
Hz due to the Bloch-Seigert effect . 



9 Experimental Results and Discussion 

The error correcting code shown in Fig. || was applied to the I^E^Ei]_, lyE^Ei^ and I^E^E']_ 
states. One can regard this code as a modification of the classical majority logic code for 
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protecting against bit flip errors about the x-axis. It follows that the I^E^E"], should not be 
affected by the gradient induced decoherence, whereas both ijE^Ei^ and I^E^E"^ will decay 
exponentially at the same rate in the absence of error correction. These exponential decay 
rates were measured by applying the gradient diffusion procedure directly to the ly and 
states, respectively. 

In Figs. and ^, the decay of the magnetization of the data spin is plotted with and 
without the error correction procedure starting from the I^E^Eij_ and IyE^Eij_ states, re- 
spectively. In both these experiments the initial slope of the error corrected curve tends to 
zero as the decoherence time t — > 0, in accord with theoretical predictions. This was quan- 
tified by using the rate 1/r of the uncorrected decay, estimated from a linear least-squares 
fit to its logarithm, to predict the corrected decay curve from equation for totally correlated 
decoherence (9exp(— t/r) — exp(— 9t/r))/8 (see Eq. (^6])). The correlation coefficient be- 
tween the measured and predicted amplitudes were 0.9845 and 0.9704 for the I^E^E'^ and 
IjE^Ei]_ experiments, respectively. In the figures, the error- corrected data points have been 
scaled so that their mean-square value is the same as that of the corresponding predicted 
points; no other free parameters were needed for these fits. 

The scatter seen in the data points, typically about ±1% of the peak intensity after aver- 
aging of 16 repetitions, can be attributed primarily to RF field inhomogeneity, particularly 
during the long period of induced decoherence. In Fig. |^, we also show the amplitudes of 
the peak starting from the I^E^Ei], data spin state, from which it is evident that it is not 
significantly affected by the overall error correction procedure. The apparent decay rate of 
^ 0.2 sec^^ can likewise be attributed primarily to RF field inhomogeneity, which produced 
a small transverse magnetization after encoding. 

Figure ^ shows the results of a similar set of experiments on the I^E^Ejj, state, with inde- 
pendent diffusion intervals for each of the three spins to implement uncorrelated decoherence, 
as previously described. As in the correlated experiments, the decay rate 1/r obtained from 
a least-squares fit to the logarithm of the decay of the I^E^E"j_ state due to this decoherence 
procedure was used to calculate the theoretical curve (3 exp(— t/r) — exp(— 3t/r))/2 (see Eq. 
(P^D), and scaled the corresponding error-corrected data to have the same root-mean-square 
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as the theoretical curve sampled at the corresponding time points. Although the experi- 
mental difficulties of juxtaposing three independent gradient-decoherence sequences without 
picking up unwanted phase shifts induces appreciable nonrandom scatter in the data points, 
the fit is once again in accord with theoretical predictions. Also shown with a dashed line is 
the decay curve expected for totally correlated decoherence, showing that these experiments 
are indeed capable of distinguishing these two cases. 

Together, these results provide strong support not only for the theory of quantum error 
correction, but also for the ability of pseudo-pure states to reproduce the dynamics of true 
pure states. This ability is important because liquid-state NMR spectroscopy provides a de- 
gree of coherent control that is presently more difficult to obtain in other quantum systems of 



comparable complexity. Although limited to about ten spins [[16], |T8[ and hence within reach 
of classical computer simulations, such an experimentally accessible paradigm for quantum 
information processing should be quite useful particularly in the study of decoherence 
It forces one to consider, and enables one to experimentally determine, actual relaxation 
superoperators, rather than working from idealized theoretical models of decoherence. The 
development of error correcting codes that can handle such real-life decoherence is a signif- 
icant challenge whose solution is likely to be applicable to quantum information processing 
in other systems. 

Conversely, the theory of quantum error correction promises to lead to new methods 
for studying molecular dynamics through the relaxation of multiple quantum coherences 



!9|, ^ ^ ^ . This may be seen, for example, by considering the second derivative of 
the error-corrected curves at t = in the uncorrelated and totally correlated cases, which are 
— and — respectively, while the corresponding inflection points occur at ln(3)r/2 
and ln(3)r/4. Moreover, since the covariances of the random fields at the three spins occur 
in our formulae for the decay of the error-corrected states, it would appear possible to 
derive these covariances directly from sufficiently detailed measurements. One of the goals 
of on-going studies of quantum error correction by the present authors is to develop codes 
which can more clearly reveal the nature of the underlying spin and molecular dynamics, for 
example by correcting for dipole-dipole relaxation. Work in these directions is currently in 
progress. 
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Appendix: Error Correction in Mixed States 

This appendix shows that error correction will not be able to correct the ffist-order decay of 
the data spin in any diagonal mixed state for the ancillae save for their pseudo-pure ground 
state, unless the state of the data spin is taken into account. Assuming that the ancillae are 
initially uncorrelated with the data spin, such a state may be written as 

p = p\ {fi++ElEl + fi+.ElEl + i^.+ElEl + ij__E''_El_) (66) 

with nonnegative coefficients satisfying + + = 1. For each term p^E^2E^3 

obtained upon expansion, the only change in the preceding results (Eq. (|43|)) is that the signs 
of the coefficients F"^ Sz are given by & , so the partial trace after error correction is 

pI = 1 + iy.2ll + (z/y2i; + u,2ll) e'{t) , (67) 

where the partial trace is initially 

p\ = I + z/,2I^ + z/y2i; + z/,2li (68) 

and now Q'{t) = 

i {F' + (/i++ + - - /i__)F2 + (/X++ - + - fi__)F^ (69) 
- (/X++ - - + f,__)F'F'F' F^23) _ 
The derivative at t = is O'(0) = 

i ((/i++ - - /x+_(c" + 2c'') - /i_+(cii + 2c'') + fi..{c'' + 2c'' + 2c'')) , (70) 

and for c^^ 7^ the only nonnegative solution to this equation together with the normalization 
condition is = 1 and = = fj, = 0. This proves that, in any such state, the 
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three-bit error correcting code (Fig. |T]) will inhibit decoherence to first order only if the 
ancillae are in a psuedo-pure state relative to the data spin. 

More generally, consider a mixed state in which the ancillae are diagonal but (classically) 
correlated with the data spin, i.e. 

p = V fi^plE%E% (71) 
= p^+E^E^ + p^_E^Ei + p^+E^ E^ + p\_ElEi , 

where = T.{m\ei=ei^=i} l^rnpln^ etc, and Xlm/^rn = 1 (/^m > 0). Letting t;,2,3 = {pl2^32ll) 
and Ce2e3 = (Pe^eB^Iz) t>e the y and z components of the data spin's constituent density 
matrices (e^,e^ G {±1}), one can show that in this case the matrix derivative at t = 
vanishes if and only if 

0++i;++ + e+_?j+_ + e_+t;_+ + e__t;__ = (72) 

where the Q^2^3 are obtained by changing the signs of and in our expression for 
= 0++. Since 0++(O) = 0, this is a system of two linear equations in the three unknowns 
—2F^ = c^^, — 2F^ = c^^ and —2F^ = c^^, which is generally solvable. Nevertheless, 
the coefficients in these equations depend on the state of the data to be protected and its 
correlations with the ancillae, which makes it impossible to use the three-bit code to protect 
unknown data in this case. 
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Figure Captions 

1. Diagram giving an overview of the error correction procedure. The three spins correspond 



to the three horizontal lines. Vertical lines connecting them are quantum gates (unitary 



operations), where "x" indicates a target spin and a control spin. The first dashed box 



is the encoding step, which consists of two successive c-JNOT's (see text). The filled box with 



the lightning bolt indicates decoherence through random rotations about the x-axis. 'I'he 



second dashed box is the decoding step, while the third is the Toffoli gate which corrects for 



decoherence to first order. The final filled triangle and square indicates that the ancillae are 



decoupled while the data spin is observed? 



|2. RF and gradient pulse sequence diagrams to implement both the totally correlated (a 



and uncorrelated (b) decoherence models. In the figure, black boxes indicate (vr)^"' RF pulses" 



and grey boxes indicate (vr)^^ pulses, while the graph on the G'^ line indicates the gradient 



polarity. The remaining symbols are defined in the main text. 



|3. Plots showing the experimental results for quantum error correction applied to the 



I^EIlEj, state (cf. Figure The "decoherence time" is the time allowed for diffusion in 
the gradient-diffusion procedure (see text). The amplitudes are of the peak due to the data 
spin, following a 7r/2 readout pulse, relative to its amplitude with no decoherence. The 
"x" symbols mark the amplitudes of the peak following the error correction procedure at 32 
equally spaced decoherence times. The "+" symbols mark the amplitudes of the peak as a 
function of the decoherence time when the gradient- diffusion procedure was applied to the 
state with no error correction at the same 32 time points. These measurements were averaged 
over 16 repetitions of each experiment (without any phase cycling). The rate of decay of the 
peak 1/r = 2.5677 sec~^ due to decoherence was obtained by a linear least-squares fit to the 
logarithm of the amplitudes of the peak from the state (lower curve, correlation coefficient 
—0.9987), after which the decay with error correction was obtained from the theoretical 
relation (9exp(— t/r) — exp(— 9t/r))/8 (upper curve, correlation coefficient 0.9873). 

4. Plots showing the experimental results for quantum error correction applied to the 



Ki^H^l state. As in 1^'igure 3|, the "x" symbols mark the amplitudes of the peak due to the 
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data spin at 32 equally spaced decoherence times (note that the amplitudes at t = 0.0145 
and t = 0.1265 sec were treated as outlyers and omitted), averaged over 16 repetitions of 
each experiment. The "+" symbols mark the amplitudes of the peak from the ly state with 
no error correction, averaged over 8 repetitions of each experiment. The rate of decay of the 
peak 1/r = 2.4200 sec~^ due to decoherence was obtained by a linear least-squares fit to the 
logarithm of the amplitudes of the peak from the ly state (lower curve, correlation coefficient 
—0.9962), after which the decay with error correction was predicted from the theoretical 
relation (9exp(— t/r) — exp(— 9t/r))/8 (upper curve, correlation coefficient 0.9867). 

5. Plots showing the experimental results for quantum error correction applied to the 
lyE^Ej. state. As in I'lgure 3| , the "x" symbols mark the amplitudes of the peak due to 
the data spin at 16 equally spaced decoherence times, averaged over 32 repetitions of each 
experiment. The decay rate of the peak 1/r = 0.2084 due to decoherence was obtained by a 
linear least-squares fit to the logarithm of the amplitudes of the peak (solid line, correlation 
coefficient —0.6022), indicating that it is not decohered by gradient-diffusion nor otherwise 
affected by the error correction procedure. 





Plots showing the experimental results for quantum error correction applied to the 




lE^ 


El state, with independent gradient-induced decoherence for each spin. As in Figure| 



[3|, the "x" symbols mark the amplitudes of the peak due to the data spin at 32 equally 
spaced decoherence times averaged over 16 repetitions of each experiment. The symbols 
mark the amplitudes of the peak from the I^E^E'^ state without encoding, decoding or 
error correction, averaged over 16 repetitions of each experiment. The rate of decay of 
the peak 1/r = 3.2931 sec~^ due to decoherence was obtained by a linear least-squares 
fit to the logarithm of the amplitudes of the peak from the data spin as a function of 
decoherence time (lower curve, correlation coefficient —0.9970), after which the decay with 
error correction was predicted from the theoretical relation (3exp(— t/r) — exp(— 3t/r))/2 
(upper curve, correlation coefficient 0.9546). The decay curve predicted for totally correlated 
decoherence as in Figure ^ is shown with a dashed line for comparison. 
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